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ABSTRACT 


In recent years, wind energy has become one of the most promising 
renewable energy sources. Various wind turbine concepts with different 
generator topologies have been developed to convert this abundant energy 
into electric power. The doubly-fed induction generator (DFIG) is currently 
the most common type of generator used in wind farms. Usually the DFIG 


generator is a wound rotor induction machine, where the stator circuit is 
directly connected to grid while the rotor’s winding is connected to the grid 
Keywords: via a three-phase converter. This paper describes an approach for the 
independent control of the active and reactive powers of the variable-speed 


DFIG DFIG. The simulation model including a 1.5 MW-DFIG driven by a wind 
MPPT turbine, a PWM back-to-back inverter and the proposed control strategy are 
PWM developed and implemented using MATLAB/Simulink/SimPowerSystems 
SimPower Systems environment. 
Simulink E l l 
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1. INTRODUCTION 

In the last decade, interest in wind energy use has grown considerably. In Europe 30-40% of newly 
installed renewable energy capacity was from wind power [1]. The U.S and China are currently the world 
leaders and dominate the global installed wind energy [2, 3]. In 2010, the world’s generation capacity from 
wind energy was 196.630 GW and reached 240 GW by the end of 2011 [4]. 

With the expansion of this renewable energy resource and its increased penetration into the 
electrical grids, wind turbine (WT) technology is currently one of the world’s fastest growing, cost-effective 
renewable energy technologies in the market [5]. 

Wind power generation is subject to fluctuations due to the intermittent nature of wind energy and 
wind speed variations. This is more evident when multiple generators are connected to a weak grid [6]. 

Many topologies for wind turbines have been designed to reduce the fluctuations of the output 
power. Doubly-fed induction generators (DFIG) are the most widely used types of generators in wind energy 
conversion systems. This topology can offset its output power to stabilize fluctuations by a factor of typically 
up to + 30%. However, this device is still small considering the range of variation in practice of the wind 
speed. Researchers have proposed that energy storage systems are a desirable choice to further mitigate the 
effects of short-term wind fluctuations [5, 6]. 

Fixed-speed induction generators are constrained to operate near the synchronous speed, because the 
frequency is imposed by the network; the rotor speed is almost constant. Variable-speed wind turbine 
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systems, on the other hand, are able to operate over a wide range of wind speeds. Furthermore, using an 
appropriate maximum power point tracking (MPPT) strategy, variable-speed wind turbines can produce 
maximum power at varying wind speed conditions. For variable-speed wind turbine systems with limited 
speed range, e.g. + 30% of synchronous speed, the DFIG can be the solution of choice [7]. 

Unlike conventional wind turbine-driven synchronous generators, the output characteristic of the 
variable-speed wind energy conversion systems (WECS) depends not only on the dynamics of the generator 
but also on the converter control strategy employed. Control strategies of DFIG have been extensively 
discussed by several authors [8, 9]. In [10], various control mechanisms for the DFIG have been reported 
using stator flux orientation in the reference frame a-p. In [11] and [12], the stator flux oriented control of the 
generator DFIG has been studied but without any control of the grid side and DC link voltages. In [13], the 
authors studied the effect of active and reactive powers variations by a decoupling control method for the 
active and reactive powers. However, the influence of different load types on the grid and power quality 
according to the loads connected to the generator DFIG has not been addressed. This paper presents a 
comprehensive simulation study of the variations of reactive power on the grid under different load 
conditions (resistive, inductive and capacitive) and analyses the influence of these loads on the DFIG power 
factor which reflects the quality of the power generated. 

The main advantage of using DFIGs in wind turbines is that they connect directly to the network via 
the three-phase stator windings and do not require additional converters. This configuration has become now 
very popular for variable-speed wind turbines [14]. This is mainly because the power electronic converter 
must handle only a fraction of 20% to 30% of the total power generated by the DFIG [14, 15]. Therefore, the 
losses in the power electronic converter can be reduced. Today, the DFIG is the most commonly used 
variable-speed machine in production units above 1 MW. 

This paper presents a study of a DFIG dual-powered generator connected with industrial loads to 
learn the impact on the quality of voltage and current by applying inductive or capacitive loads. 

The aim of this paper is to design an indirect vector control decoupling strategy to achieve 
independent control of the active and reactive powers of the DFIG-WT conversion system. A detailed model 
of the DFIG is used to design and evaluate the proposed control strategies. The authors present a technique of 
the stator quadrature flux along the axis to control the active and reactive powers of the DFIG. The proposed 
control approach is evaluated under resistive, inductive and capacitive loads connected between the DFIG 
and the grid. 

The rest of the paper is organized as follows: Section 2 introduces the overall simulation model 
structure including mathematical models of the wind turbine. The proposed MPPT method used in the model 
is also presented. Section 3 of the paper presents the model of the DFIG and its indirect vector control 
scheme using Park transformation, the control of the DC side voltage. Finally, the simulation results and 
conclusions are summarized in Section 4 and 5 respectively. 


2. PROPOSED WIND ENERGY CONVERSION SYSTEM MODEL 

The proposed circuit is shown in Figurel. It consists of a three-blade wind turbine connected to a 
variable-speed DFIG. The stator circuit of the DFIG is directly connected to the grid, whereas the rotor 
winding is connected to the grid via two PWM-controlled IGBT-converters. The two converters are coupled 
through a DC link capacitor. 

The rotor-side converter controls simultaneously the active and reactive powers by adjusting the 
amplitude, frequency and phase of the rotor voltages. The aim of the grid-side converter is to regulate the DC 
link voltage. The inputs to the control system are the voltages and currents on the source side and the voltage 
on the DC side. These quantities are transformed into their d and q components. A phase-locked loop circuit 
is used to synchronize the frequency of the system with the network frequency. 
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Figurel. Wind turbine with a DFIG 


2.1 Modeling the wind turbine with MPPT strategy 


The relationship between the wind speed and the aerodynamic mechanical power extracted from the 
wind can be described as follows [16, 17]: 


1 
P= ZOAR L, (A, B).v 


m 


(1) 


Pm: is the mechanical power of the wind turbine [W], 

B: is the orientation angle of the blades [°]. 

The power coefficient Cp defines the aerodynamic efficiency of the wind turbine. It depends on the 
characteristic of the turbine and is a function of the speed ratio à and the orientation angle B of the blade: 








C 2s 
C,(A, B)=C, (S-(c,)-C | “+C,A 
(2) 
The values used for the coefficients C1 - C6 are given in Table 1 [18]: 
Table 1 Coefficients of the turbine 
C, C, C; C4 C; Co 
0.5176 116 0.4 5 21 0.0068 
With 
d. 7 l E 0.035 
A A+0.08f8 B +l (3) 
À is defined as the ratio of the linear velocity at the end of the blade and is given by: 
a Pak 
k (4) 
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It can be observed from Figure 2 that the power coefficient reaches a maximum for a pitch angle of 
0° and a particular value Anom = 8.1 of the velocity ratio. The power coefficient value corresponding to 
Anom is Cpmax = 0.48. 


C 0B) 








Figure 2. Characteristic of the power coefficient as a function of à 


Depending on the wind aerodynamic conditions, there exists an optimal operating point which 
allows the maximum power to be extracted from the turbine. This can be achieved by either controlling the 
rotational speed of the turbine or the power of the turbine. Several MPPT methods can be used either with or 
without the knowledge of the wind turbine characteristics [19]. In this study, the first method has been 
applied. The optimum rotation Qm, opt for the mechanical transmission of the maximum wind turbine is 
given by [20]: 





m, opt R (5) 


the following relation can be deduced: 


_ 3 
P max K p op Ym, opt (6) 
where: 
1 R 
K pon = 5 PC 43 
,opt 9) p,max A’ 
bi (7) 


thus, the corresponding optimum torque is: 


T — re — K 2 
m,opt Q —— => p, opt “m, opt 


m, opt (8) 





On the power characteristic of a turbine (Figure 3), the locus of the point representing the maximum 
power is obtained by adapting the speed of the turbine (thick curve) to that of the wind speed. 
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Figure 3. The characteristics of system optimum turbine 


2.2 Modeling of the DFIG 
The Park model of the DFIG is given by the following set of equations [21]: 











. do 
Vas = R las TO, qs + dt 
. do 
Vis = Ris + 0.9,, + 7 
do 
v, =Ri,-(@,-@ + 
dr r-dr ( S r )¢,, dt 
d 
v,, = Ri, +(@,-@,)¢,, + Par 
dt (9) 
The stator flux equations are: 
a T Lias F Lpnlar 
És = Lis + Libor (11) 
Similarly, the rotor flux equations are: 
[r = Lig T Li las 
Pip = Lyle + Labgs (12) 


In these equations, Rs, Rr, Ls and Lr denote respectively the resistances and inductances of the 
stator windings and rotor, Lm is the cyclic mutual inductance, wr = P.Or is the rotor speed (with P the 
number of pole pairs) and ws is the synchronous angular speed. vds, vqs, vdr, vqr, ids, iqs, idr, iqr, @ds, øqs, 
édr and gar are respectively the direct and quadratic components of voltages, currents and fluxes in the stator 
and rotor. The active and reactive powers of the stator and rotor are obtained as: 


| va tg Tra das ) 
Q, = iv, las — Vas lgs ) 


(13) 
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Li =(v, Lay, + Vay Lay ) 


3 


0 = 2 (va Ly, — Var La, ) (14) 


The mechanical and electromagnetic torque are given by the following equations: 


js Se ne 
dt (15) 


| aia = (¢,, Li, + Pi, Ly, ) 
; (16) 


Where J is the moment of inertia and f is the viscous friction coefficient. 


3. CONTROL SCHEME OF THE GENERATOR AND CONVERTERS 
3.1. Control of the Rotor Side Converter 


The stator flux vector is chosen to be aligned with the q-axis Park reference frame. The grid is 
assumed to be stable and therefore øsq is constant. The resistance Rs of the DFIG stator is neglected. 








Pis = 0 
V S 
Ds m D, oo P 
s (17) 
5 = V, 
v, =O 
4s (18) 
Where Vs represents the r.m.s. value of the grid voltage. The torque (16) becomes: 
PL 
L Z -5 slar 
B (19) 
The stator flux equations (11) become: 
p A Liis nG Lidar 
=Li +L i 
Pys slqs im! gr (20) 
From (20), the equations linking the stator and rotor currents are: 
ee 
las =~ L lar 
s Ln 
qs L L qr 
S S (21) 


The active and reactive powers of the DFIG given by (18) are re-written as: 
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3 3 L 
be a 
2 2L, 
3 3/ VL V? 
Q, = -Vi S = 2 = L, aa > 
2°" 2 L ” aL 
S S S (22) 
Then the rotor voltages can be written as: 
| L, 
Var = Yi 7 Olip Qs L ii L,o lor 
Vor — Vor T Olip L,o Lay (23) 
| di 
v, =Ri, +Lo—& 
dr r dr r dt 
! ° di r 
Vv, =R,i,, + L,o— 
dt (24) 


Where the slip angular velocity slip, and total leakage coefficient o are given by: 


2 


m 





O stip =0, TO,» a 


S r 


The rotor side controller consisting of the active and reactive power controllers is shown in Figure 4. 


œslip(Lm/Ls.øqs+0oLridr) 






Current 


(Fast Control 
Loop 






œslip o Lr 
idr 


idqr kA trabc 4 


oslip oslip 


encoder 


= Or 
THS 
+ 


“a © 


Poga Calculate active vdqs 
power based on dq idqs isabc A 
Omeds components 





Transformer 


Three phase 
Grid 


Figure 4. Rotor side controller for the DFIG 


3.2 Control of the grid side converter 
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The grid side converter can also be used to control the power factor of the system by adding a 
control loop for the reactive power via the q-axis current. In the system, the quadratic component of the 
reference current is set to zero in order maintain a unity power factor. A similar control strategy was applied 
to the grid side converter, since the active and reactive powers of the grid side converter can be controlled 
independently by acting on the d and q components of the grid side voltage [5]. The grid side converter 
model in Park coordinates system is [22-25]: 


{ — 0 : 
H L; He l i = 
dt l, ae R; l L; VG V 


a (25) 


Where: Rf and Lf are respectively the resistance and the leakage inductance of the grid-side transformer, vd 


and vq are the source voltage components and voq vod are the inverter voltages. Using the decoupling 
method by compensation, the inverter voltages can be written as: 


Yva Bae Vai 
Vog = Eog “Yogi 


k (26) 
With the control coefficients: 
di 
Koai — /¢ —< 
dt 
di 
v,,=L,— 
oql f 
dt (27) 
And the coefficients of compensation are obtained as: 
a 
Cog = -L Ol, E Ri p Va (28) 


Neglecting the converter losses, the DC bus voltage vdc varies with the power exchanged between the 
turbine and the network and is given by: 


dv, -r 


=] =f l 
dc dc _ rotor dc _ grille 
dt 


C 
(29) 


The proposed control strategy for the DC link is depicted in Figure 5. It includes two control loops, 
an inner loop and an outer loop. The outer loop consists of a proportional and integral (PI) controller for 
regulating the DC link voltage. The output of the DC voltage controller represents the d-axis current 
reference idref from the source which is then compared with the measured current Idmeas. The quadrature 
component of the source current iq is used to control the flow of reactive power. As discussed earlier, the 
reactive power reference is set to zero in order to obtain a unity power factor. 
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Figure 5. Controller side of the grid DFIG 


4. RESULTS AND DISCUSSION 

Figure 6 shows the overall model of the WECS and control scheme implementation in 
MATLAB/Simulink/SimPowerSystems. In the following, some simulation tests are presented to illustrate the 
performance of the proposed control system for the WECS. 





Control System DFIG Detailed Model | | es 


Figure 6. MATLAB/Simulink diagram of DFIG and control scheme 
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The Side Grid Converter (SGC) adjusts the amplitude and frequency of the signal to be sent to the 
rotor of the DFIG: It acts on the voltages at the terminals of the rotor circuits. It varies the speed of the 
turbine and thus the power extracted. 

Figure 7 shows the responses of direct (ird) and quadrature (irq) DFIG rotor currents and the stator 
active (Ps) and reactive (Qs) powers. The direct component ird controls the active power Ps. For a nominal 
power of Ps = 1.5 MW, the required forward current reaches a value of 1.5 kA. The reactive power Qs is 
controlled by the reactive component of current irq. The results show a good tracking performance when the 
reactive power reference is stepped from 0 to O.SMVAR and then -O.5MVAR. 

The reactive power supplied to the grid can be controlled by the reactive power generated or 
absorbed by the Rotor Side Converter (RSC) connected to the rotor. The reactive power is exchanged 
between this converter and the network, through the generator. Indeed, it absorbs reactive power to 
compensate for mutual inductances and leakage inductances. The converter connected to the network can 
also operate as a reactive power compensator. 

It is noted that any variation in the rotor or rotor voltage frequency has a direct influence on the 
power and torque. 

The quadrature component of the rotor current controls the active power, and the direct component 
controls the reactive power exchanged between the stator and the network as shown in Figure 7. It can be 
observed that both active and reactive powers of the DFIG follow their references. 

The DC link voltage, the amplitude modulation index (MI) and the voltage vra and current ira in 
phase A of the rotor are shown in Figure 8. It can be noted that the DC voltage tracks perfectly its reference 
of 1.5 kV. 

Figure 9 shows the responses of the electromagnetic torque, rotor angular speed, the speed ratio À 
and the power coefficient Cp. The power factor Cp is expected to reach an optimal value of 0.48 after a short 
transient while the speed à reaches only the maximum of 8.1. 
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The DFIG is now connected to resistive, inductive and capacitive loads respectively. The simulation 
results shown below are obtained at the same wind speed of 11 m/sec. Figure 10 shows the responses of the 
real (Ps) and reactive (Qs) powers of the stator, the real (PL) and reactive (QL) powers of the load and the 
resulting currents ird and irq in the case of a resistive load. Note that the grid absorbs the DFIG 1.5 MW 
active before the application of the load. After application of a load of 1.5 MW at time t = 1 sec, the power 
supplied by the DFIG is then consumed by the load and the power of the grid becomes zero. The results are 
illustrated in Figure 11. It can be noted that the voltage vdc follows its reference of 1.5 kV. The grid and load 
currents reached their steady-state values of 1.5 kA and O A respectively before the application of the load. 
After the application of the load, the current of the load becomes 1.5 kA and that of the grid becomes 0 A. 
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Figure 11. DC link voltage, amplitude modulation Figure 12. Active and reactive powers of stator, load 
index (MI), phase A current of the grid and load and and rotor currents (inductive load conditions) 


the current of DFIG rotor (resistive load conditions) 


Figure 12 shows the responses of real and reactive powers of stator (Ps and Qs) and those of the 
load (PL and QL) and the responses of the DFIG rotor currents ird and irq for an inductive load (PL =1.5 
MW and QL =0.5 MVAR). Note that the grid absorbs the active power generated by the DFIG before 
applying the load. After application of the load at the instant t = 1 sec, the necessary power to be consumed 
by the load is supplied by the DFIG and the power of the grid becomes zero. The load, the DC link voltage 
vdc, modulation index (IM), the currents of the grid and load and the DFIG rotor currents are illustrated in 
Figure 13. Similar remarks and conclusions can be drawn as in the case of resistive load. 
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Figure 13. DC link voltage, amplitude modulation Figure 14. Active and reactive powers of stator and 
index (MI), phase A current of the grid and load load and rotor currents (capacitive load conditions) 
and the current of DFIG rotor current (inductive 

load conditions) 
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A capacitive load is now introduced (PL = 1.5 MW and QL = 0.5 MVAR) and Figure 14 show the 
results. Note that, in this case, the grid absorbs the active power generated by the DFIG before the load is 
applied. After application of the load at the instant t = 1 sec, the power supplied by the DFIG is consumed by 
the load and the power of the grid becomes zero. Figure 15 shows, the DC link voltage vdc, the modulation 
index (IM), the currents of the grid and the load and DFIG rotor currents. Table 2 summarizes the results 


obtained with the three types of loads and the effects on the power factor of the DFIG. 
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Figure 15. DC link voltage, amplitude modulation index (MI), phase A current of the grid and load and the 


current of DFIG rotor current (capacitive load conditions) 


Table 2 Comparison of the three types of loads. 


load active power of DFIG [MW] reactive power of DFIG [MVAR] Power factor CosQpric 
resistive 1.5 0 1 

inductive 1.5 -0.5 0.95 lagging 
capacitive 1.5 0.5 0.95 leading 


Figure 16 shows the current waveforms of phase A for the three different types of load and the grid 
phase A voltage. This figure clearly illustrates that the voltage and the current generated by the DFIG are in- 
phase when the load is resistive. The current is lagging and leading the voltage when the load is respectively 


inductive and capacitive. 


One can say that the shapes of the curves show an acceptable tracking of the stator reactive power. 
The same applies to the stator active power and the total reactive power supplied to the network. This method 
of control presents remarkable performance in terms of reactive power control in the network through the 
wind system. This setting allows us to have a negative reactive power (capacitive behavior) or positive 


(inductive behavior) and even zero (unity power factor). 


It can be noted that the reactive power that can be generated or absorbed is limited by the level of 
the rotor currents, these currents are also imposed by the switches of the IGBT power devices. 
Figure 17 shows the waveforms of the instantaneous active power generated by the DFIG under 


these three load conditions. 


In the case of capacitive load, ripples can be observed in Figure 17 due to the stator voltages of the 
DFIG. On the hand, in the case of inductive load, voltage drops appear in the network voltages which caused 
a reduction of these ripples in the stator power waveform. In addition, the inductance can filter the line 


currents which leads to a significant improvement of the power quality. 
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Figure 16. Phase A grid voltage and current Figure 17. DFIG stator active powers for resistive, 
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5. CONCLUSION 

This paper presented the modeling, simulation and control of a DFIG based wind conversion system 
connected to the grid. The overall system was implemented using MATLAB/Simulink and 
SimPowerSystems toolbox. The indirect vector control method is used to control the stator flux and hence the 
active and reactive powers of the DFIG. The reference of the active power generated by the DFIG is set 
according to the wind speed using a maximum power point tracking (MPPT) strategy. The reactive power 
set-point is varied according to the type load connected to the generator. Three types of load scenarios have 
been considered (resistive load inductive and capacitive) to demonstrate the influence of the load on the 
DFIG. The simulation results show the effectiveness of the decoupling controller and MPPT strategy under 
these three load conditions. The grid-connected DFIG simulation model presented in this work can be easily 
extended to include other compensating devices such as SVC (Static Var Compensator) or ASVC (Advanced 
Static Var Compensator) and can also be used to investigate other more robust control strategies to further 
improve the performance of the system in response to grid disturbances. 
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